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ABSTRACT 

We use SDSS-DR4 photometric and spectroscopic data out to redshift z 0.1 combined with 
ROSAT All Sky Survey X-ray data to produce a sample of twenty-five fossil groups (FGs), defined as 
bound systems dominated by a single, luminous elliptical galaxy with extended X-ray emission. We 
examine possible biases introduced by varying the parameters used to define the sample and the main 
pitfalls are discussed. The spatial density of FGs, estimated via the T^/Vmax test, is 2.83 xlQ~^h^^ 
Mpc~'^ for hx > 0.89 X lO^'^/i^g^ erg s~^ consistent with Vikhlinin et al. (1999), who examined an 
X-ray over luminous elliptical galaxy sample (OLEG). We compare the general properties of FGs 
identified here with a sample of bright field ellipticals generated from the same dataset. These two 
samples show no differences in the distribution of neighboring faint galaxy density excess, distance 
from the red sequence in the color-magnitude diagram, and structural parameters such as a.^ and 
internal color gradients. Furthermore, examination of stellar populations shows that our twenty-five 
FGs have similar ages, metallicities, and a-enhancement as the bright field ellipticals, undermining 
the idea that these systems represent fossils of a physical mechanism that occurred at high redshift. 
Our study reveals no difference between FGs and field ellipticals, suggesting that FGs might not be a 
distinct family of true fossils, but rather the final stage of mass assembly in the Universe. 
Subject headings: Galaxies: formation - Galaxies: evolution - Galaxies: fundamental parameters 



1. INTRODUCTION 

Searching for preserved remnants of physical processes 
that occurred in the cosmic past is one of cosmologists' 
main tools in developing our understanding of how galax- 
ies, groups, and clusters were formed and evolved, and 
ultimately came to define the large scale structure we 
observe. Recent studies have shown that there seems 
to be a class of systems with exceptional preservation, 
known as fossil groups (FGs, Ponman et al. 1994). These 
systems consist of isolated, luminous early- type galaxies 
embedded in an extended X-ray halo. The main moti- 
vation for seeking such galactic systems stems from the 
fact that the merger time scales for L > L* group galax- 
ies is much shorter than the cooling time scales for the 
hot gas component within which these bright galaxies are 
embedded (e.g. Barnes 1989, Ponman & Bertram 1993). 
In this simplistic view FGs are a natural consequence of 
the merging process in normal groups and therefore can 
be used to trace the history of coalescence in the not-so- 
high-redshift Universe. 

From the observational viewpoint, FGs are seen in the 
optical as large elliptical galaxies, but with X-ray lumi- 
nosities comparable to those of an entire group of galaxies 
{Lx > 4.4 X lO^^ergs/s; see Jones et al. 2003, here- 
after JO03). To date, no clear picture of their origin 
has emerged from the collected data. Two main scenar- 
ios have been suggested: FGs result from the complete 
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merging of galaxies that constituted a loose group in the 
past, collapsing at an early epoch, but never becoming 
incorporated into clusters (Hausman & Ostriker 1978; 
Ponman et al. 1994; Jones et al. 2000; Khosroshahi, 
Jones & Ponman 2004); or they originate from a region 
that inhibits formation of L* galaxies in these groups 
leading to an atypical galaxy luminosity function (LF, 
Mulchaey & Zabludoff 1999). 

These systems have been detected out to redshifts of at 
least z ~ 0.6 (Ulmer et al. 2005), have high M/L ratios 
suggesting low star formation efficiency (Vikhlinin et al. 
1999) and represent 8% to 20% of all systems with simi- 
lar X-ray luminosities. Ten FGs are already known and 
well studied from optical and X-ray data (Mulchaey & 
Zabludoff 1999, Romer et al. 2000, hereafter ROO, Jones 
et al 2003, Khosroshahi, Jones, & Ponman 2004, Sun et 
al. 2004, Yoshioka et al. 2004, Ulmer et al. 2005). A 
recent paper by Santos et al. (2007) presents a sample 
of 34 FG candidates based on SDSS data. In this work, 
we use optical (SDSS) and X-ray (ROSAT All Sky Sur- 
vey - RASS) data to define a sample of FG candidates. 
We select FGs following a strategy similar to Jones et 
al. (2003), using spectroscopic (age, metallicity, a/Fe 
enhancement) and photometric (a4, color gradient) pa- 
rameters to further constrain what should be a FG. The 
main goal is to establish a link between FGs, compact 
groups and isolated ellipticals based on their global prop- 
erties and simple expectations, if mergers are the domi- 
nant events determining the evolution of a FG. 

This paper is organized as follows. Section 2 presents 
the galaxy catalog used to search for FGs, describing the 
parameters and the X-ray measurements that define the 
search. The actual selection of FGs is then detailed. Sec- 
tion 3 presents the control sample used to compare FG 
properties to those of "normal" ellipticals, which in prin- 
ciple did not result from merging. Section 4 discusses 
the magnitude distribution of the brightest and second 
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brightest galaxies of FGs, while in Section 5 we charac- 
terize the distribution of the faint galaxies around FGs. 
Since FGs are thought to be old systems, in Section 6 we 
examine the colors of FGs relative to the red sequence 
of early type galaxies. Section 7 deals with the inter- 
nal structure of seed galaxies of FGs, as measured by 
their structural parameters and internal color gradients. 
In Section 8, we analyze stellar population properties of 
FGs, by looking at the distributions of age, metallicity, 
and a enhancement. Finally, in Section 9, we summa- 
rize and discuss the main results. Throughout the pa- 
per we use a cosmology with fl^ = 0.3, Qa — 0.7, and 
Ho = 75kms-iMpc-^ 

2. SELECTION OF FGS 

FGs are defined as bound systems of galaxies asso- 
ciated with extended X-ray sources. The total optical 
luminosity of these systems is dominated by a bright el- 
liptical galaxy, and their LF exhibits a gap between th e 
first and second rank galaxies (see e.g. lJones et al.l[2003f ). 
The sample of FG candidates defined here is selected on 
the basis of spectroscopic and photometric data from the 
Data Release 4 (DR4) of the Sloan Digital Sky Survey 
(SDSS), and on the basis of RASS X-ray imaging. The 
optical selection is described in Sec. 12. 1[ while the X-ray 
measurements and the definition of the FG sample are 
presented in Sec. 12.21 

2.1. The SDSS catalog and the optical selection 
parameters 

The optical selection is performed on a complete 
volume-limited catalog of 91563 galaxies, retrieved from 
the SDSS-DR4 database ^ through the CasJobs facil- 
ity The catalog consists of all galaxies with abso- 
lute r-band magnitudes Mr brighter than —20 and with 
spectroscopic redshifts between 0.05 and 0.095. Abso- 
lute magnitudes are obtained from the SDSS r-band Pet- 
rosian magnitudes, corrected for galactic extinction using 
the re ddening maps of Schlcgcl, Finkbcincr and Davis 
()1998f l. The lower redshift limit o f the catalog is cho - 
sen to minimize the aperture bias (jGomez eTaDHool), 
which strongly affects large nearby galaxies, while the 
upper redshift limit guarantees a high level of com- 
plet eness, estimated thro ugh Schmidt's V/Vmax test 
(see lSorrentino et al.ll2006l ). The magnitude limit of —20 
corresponds approximately to the apparent magnitude 
limit of the SDSS spectroscopy at redshift z ^ 0.095 
(r ~ 17.8). 

Since FGs are systems characterized by a gap in the 
galaxy LF, we first select optical candidates by search- 
ing the SDSS spectroscopic database for pairs of galax- 
ies whose magnitude difference is larger than a given 
value, AMmin- For each galaxy in the catalog, we se- 
lect all objects whose projected distance on the sky from 
the target galaxy is smaller than a maximum radius, 
-Dmax, and whose redshift difference is smaller than a 
given value, A(c2:max)- In other words, for each galaxy, 
we select companion galaxies in a cylinder centered on 
the galaxy, with a radius of I?max and a semi-height of 
A(czmax) along the line of sight. The value of i^max is 
fixed in physical units (Mpc), and the selection is done 

http://www.sdss.orK/DR4l 

* http://casjobs.sdss.Org/casjobs/l 



by transforming the value of -Dmax to the apparent size 
that corresponds to the redshift of the seed galaxy. A 
bright galaxy is defined as a possible FG if (i) there is at 
least one companion selected through the above proce- 
dure, and (ii) if all the companions are fainter, with the 
magnitude difference greater than AMmin. Thus, our 
optical FG candidates are selected on the basis of three 
parameters: A(czi„ax), -Dmax, and AMmin. The pa- 
rameters A(czmax) and i^max are used to select the 
companions of a given galaxy in the SDSS catalog, while 
the parameter AMmin characterizes the gap in the LF 
between the galaxy and its companions. 

To determine suitable values of A(c2;max), -D^ax 
and AMmin, we examine how the number of optical 
FG candidates and the corresponding contamination rate 
change as a function of these parameters. The param- 
eters are chosen in such a way that the contamination 
rate is at most 15%, defined as the fraction of FG candi- 
dates that would be selected (for a given set of values of 
A(cz,„ax), -Dmax and AAfmin) from a random distribu- 
tion of galaxies with the same mean galaxy surface den- 
sity, magnitude, and redshift distributions as the SDSS 
catalog. We generated randomized galaxy c atalogs by 
applyi ng the shuffling method, as described in iGal et all 
()2003l ). First, we bin the SDSS catalog in 4° in right as- 
cension and declination and Az — 0.005 in redshift. The 
randomized catalog is then obtained by shuffling galax- 
ies in each bin by randomly exchanging their right as- 
cension, declination and redshift. The above bin sizes 
are chosen such that both the mean surface density and 
the redshift distribution of galaxies in the SDSS cata- 
log are preserved, without overly reducing the number of 
galaxies in each bin. We verified that varying these bin 
sizes by 20% does not change the contamination rate es- 
timates. Fig.[T]compares the characteristics of the galaxy 
catalog with those of one of the shuffled catalogs. The 
shuffled catalogs reproduce well the sky coverage of the 
SDSS-DR4 catalog as well as the corresponding distribu- 
tions of galaxy magnitudes and redshifts. Since the shuf- 
fling procedure washes out small scale structures from 
the galaxy catalog, the FG candidates that are selected 
from the shuffled catalogs only arise as rare and random 
instances in the magnitude, redshift, and spatial distribu- 
tions of the galaxies, instead of a population of objects 
that follow from some physical evolutionary path. To 
reduce statistical fluctuations in the contamination rate 
estimates, we repeat the shuffling procedure 100 times 
and estimate the contamination rate as the ratio of the 
number of optical FG candidates to the mean number of 
candidates selected from the shuffled catalogs applying 
the same A(czmax), i'max, and AA/min criteria. 

Figs. and [3] show the number of FG optical 
candidates and the contamination rate as functions 

of A(cZmax), -Dmax and AMmin- We vary A(cZmax) 

from 0.00025 75 km s~^), corresponding to ~ 2.5 
times the typical redshift uncertainty of the SDSS spec- 
troscopic sample ^, to 0.0065, corresponding to ^ 2000 
km s~^, i.e., about two times the typical velocity disper- 
sion of a rich cluster of galaxies. The minimum value 
of -Dmax is chosen to be 0.3Mpc, which is about twice 
the typical core radius of galaxy clusters (jAdami et al.l 

^ For the SDSS-DR4 ma in galaxy sample, the typi cal redshift 
accuracy is 30 km s~^ (see |http:/ /www.sdss.org/dr4/| 
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Fig. 1. — Comparison of the characteristics of the galaxy cat- 
alog with those of one of the shuffled catalogs. The upper panel 
plots declination versus right ascension and compares the sky area 
covered from both catalogs. The lower left and right panels show 
the redshift and magnitude distributions, respectively, of galaxies 
in both catalogs. As shown in the upper panel, galaxies from the 
shuffled and the galaxy catalogs are plotted with grey and black 
colors, respectively. 
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Fig. 2. — Number of FG candidates as a function of the selection 
parameter Azmax- The panels correspond to different values of 
AMinin (shown in the upper right of each plot), while different line 
styles correspond to different values of Dmax, as shown in the 
lower right corner of the figure. 



at A(c2:niax)~ 0.001 for a wide range of values of both 
Dma.x and AMmin- Hence, to maximize the sample size, 
we adopt this value of A(czniax)- For Dmax, we see that 
for fixed values of A(c2;„iax) and AMmin the number 
of candidates does not change significantly for 0.35 < 
Dmax'^ O.SMpc. On the other hand, the contamination 
rate changes by a factor of 2 over the full range of values 
of Dniax- We choose a suitable value of Dmax by fixing 
an upper limit to the contamination rate. Examining 
Fig. [3] and considering all the values of AMmin that cor- 
respond to A(czmax)= 0.001, contamination rates less 
than 15% are achieved only for Dmax^ 0.35Mpc. Since 
the number of FG candidates decreases as a function of 



Dmax, we adopt D 



0.35Mpc. Finally, we have to 



determine the choice of AAfmin- Fig.|4]plots the contam- 
ination rate as a function of AMmin for all possible values 
of Dmax and A(c2;max)- We see that the contamination 
is almost constant for AMmin< 2.2, with an increase for 
larger values of AMmin- Fig. [5] plots the number of FG 
candidates as a function of AAfmin for A(c2;max)= 0.001 
and Dmax= 0.35Mpc. The number of FG candidates 
rapidly decreases as AMmin increases, reaching almost 
at AMmin~ 2.5. In order to obtain a statistically signif- 
icant number of FG candidates without overly reducing 
the AAfmin gap, we decided to adopt AMmin= 1-75 mag. 

To summarize, the values of the three parameters 
A(czmax), Dmax and AAfmin are chosen as a com- 
promise between the number of selected FG candidates 
and the corresponding contamination rates. We adopt 
the following parameters: A(czmax)= 0.001, Dmax= 
0.35Mpc and AMmin = 1.75mag, resulting in a hst of 
578 optical FG candidates. Note that the Dmax value is 
close to half of the virial radius (0.37 hyg^Mpc) assuming 
a te mperature of 1 Key and the mean redshift of the sam- 
le (lEvrard et al.iri996[ ). As found bv iKhosroshahi et al.l 
20071 ). based on Chandra data, a typical temperature 
for FGs is 1 Kev, with a few systems having a tempera- 
ture as high as ~ 3 Kev. In the latter case, half the virial 
radius is 0.64 Mpc. Using such a large Dmax results in a 
high contamination rate for whatever value of A(czmax) 
and AMmin (see Fig.Hl). 

We note that the above selection of FG systems can be 
affected by the SDSS spectroscopic incompleteness, be- 
cause of the SDSS fiber collision limit preventing neigh- 
boring fibers from being closer than 55". This prevents 
complete spectroscopic coverage of objects in dense re- 
gions, even in cases where multiple tiles overlap. It is 
therefore, possible to have galaxies adjacent to some 
optical FG candidates which would nominally disqual- 
ify them, but lacking SDSS spectroscopy, the neighbors 
would not be taken into account by our selection proce- 
dure. This issue is addressed in Sec. 12.2.21 



Il998f l. Furthermore, we consider only values of Dmax 
smaller than 0.6Mpc, since higher values yield contam- 
ination rates greater than ~ 20% for all possible val- 
ues of the other two parameters. Finally, we consider 
only values of AA/min in the range of 1.5 to 2.5 mag. 
The value of 1.5 is chosen to have a reasonable mini- 
mum gap in the LF of FG candidates, while values of 
AMmin larger than 2.5 are excluded since they overly 
reduce the number of candidates. Looking at Fig. [21 we 
see that the number of candidates shows a maximum 



2.2. X-ray Analysis and the Final Sample of FGs 
2.2.1. Measuring X-ray Fluxes with RASS 

The X-ray luminosity Lx of each optical FG candidate 
is estimated using RASS countrates^°. The first step is 
to estimate the background contribution, by randomly 
selecting one hundred boxes each with an area of 100 
square arcminutes, located throughout the 6.5° x 6.5° 
field containing the target, and then taking the median 

^'^ ftp:/ /ftp.xray.mpe.mpg.de/rosat/archivc/900000 
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Fig. 3. — Same as Fig.|2]but showing the corresponding contam- 
ination rates as estimated through the shuffling method (see the 
text). In each panel, the dotted horizontal line marks the upper 
limit for the acceptable level of contamination in the sample of FG 
candidates. 
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Fig. 4. — Contamination rate as a function of AAf,„in. Each plot 
corresponds to a different value of A(cZniax)i as shown in the upper- 
right of each panel. Curves with different line styles correspond to 
different values of Dmax, and are the same as in Fig. [2] 



and rms of these hundred values. The count rates are 
integrated in three apertures, with radii of 5, 10, and 
20 times the effective radius, re, of the FG seed galax- 
ies. Effective radii are measured from the r-band SDSS 
images, as described in App. [Cl We considered a detec- 
tion as significant when the S/N is at least 3-cr above 
the background, otherwise we derive upper limits as in 
Beuing et al. (1999). We selected for further analysis 
only those systems with a significant X-ray detection in 
at least one of the three measured apertures, resulting in 
a sample of 113 (out of 578) FG candidates. 

To convert the measured total count rate into an un- 
absorbed X-ray flux in the ROSAT 0.5 - 2.0 keV energy 
band, we use the HEASARC PIMMS tool. We assume 




Fig. 5. — Fraction of FG candidates in the SDSS catalog as 
a function of AAfmin for the values of Dniax= 0.35Mpc and 
A(czinax)= 0.001 adopted to select FG candidates. The black 
curve is obtained from a cubic interpolation of the data. 



a Raymond-Smith spectrum to represent the hot plasma 
present in the intracluster medium (Raymond & Smith 
1977) with solar metallicity and an interstellar hydrogen 
column density along the line-of-sight given by Bajaja 
et al. (2005). The plasma temperature is fixed at 1 keV, 
which is typical for very bright elliptical galaxies (Beuing 
et al. 1999), although the dependence of the flux on the 
temperature used is insignificant (< 5%). Finally, the 
rest frame 0.5-2.0 keV flux was obtained by applying the 
K-correction as in Bohringer et al. (2000). 

2.2.2. Contamination by Spirals, AGN, and Clusters 

To refine the selection of FGs, we considered four pos- 
sible sources of contamination: (i) spiral galaxies; (ii) 
AGN; (iii) superpositions with rich clusters of galaxies. 

(i) Since FGs are dominated by a bright elliptical 
galaxy, we removed contaminating spiral galaxies from 
the list of optical candidates. We fit the SDSS r-band im- 
age of each galaxy with a seeing convolved Sersic model 
(see App. [Cl for details), and flagged as spirals those ob- 
jects with a detectable spiral arm pattern in the resid- 
ual map obtained after model subtraction. This proce- 
dure eliminates ninety-one (out of 578) contaminating 
systems. 

(ii) We flagged those FG optical candidates with AGN 
spectral signatures and excluded them from the analy- 
sis . AGN were selected f ollowing the criteria described 
in ISorrentino et all (|2006D . We removed fourty-five (out 
of 578) galaxies. 

The above selections leave a sample of 102 FG can- 
didates with significant X-ray detections, no significant 
AGN signatures, and early- type morphology. 

(iii) The resolution and S/N of RASS is insufficient to 
differentiate between X- ray emission from a real FG and 
that from a rich cluster which happens to be along the 
line of sight. Moreover, for a rich cluster, the SDSS spec- 
troscopic incompleteness at the bright end might mimic 
the lack of bright companion galaxies required to be a 
FG. Thus, we rejected FG candidates within 1.5 Mpc (at 
the FG redshift) from a rich Abell cluster. We discarded 
i? > Abell clusters with the exception of Abell 690 
(richness R — 1; FG22 of Tab. 2) for which we veri- 
fied that there are no galaxies with available photometry 
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Fig. 6. — Diagnostic diagrams showing the nine FG candidates 
(filled circles) which still show some emission lines in the SDSS 
spectra. The empty circles represent the six field galaxies (FS; see 
Sec. [Sj with some residual emission features in the SDSS spectra. 
We also display the regions defining the locus for AGN (Seyfert and 
LINERs) and starbursts (HII). Notice that the figure also shows 
the transition region between AGN and HII (Comp). 



from SDSS that might invahdate our adopted fossilness 
definition. Twenty-eight (out of 102) FG candidates were 
discarded with this criterion, leaving a sample of seventy- 
four FG candidates. 

Although we eliminated AGN and spiral galaxies from 
the sample, we further investigated the presence of elu- 
sive emission lines in the SDSS spectra of the 102 FGs 
with significant X-ray detections. For each spectrum, we 
modeled the absorption lines with a combination of stel- 
lar population models using the STARLIGHT code (Cid 
Fernandes et al. 2005). Further details of the proce- 
dure are presented in Sec. [H After subtracting the old 
stellar component, we measured the remaining Oj, On, 
OiUjHa, Hf3, and Nu emission features. Fig. [5] displays 
the diagnostic diagrams defined with these indices. Re- 
gions corresponding to different types of active galaxies 
are pl otted according to the definition of iKewlev et al.l 
((200l . Only nine FG candidates have residual emis- 
sion lines typical of AGN. In agreement with expectations 
from our morphological selection, no starburst galaxies 
(HII) arc found in the FG sample. Moreover, almost all 
of the AGN are classified as LINERs, w hich have a typ- 
ical X -ray luminosity of 10"*° ergs/s (see iKomossa et all 
Il999l two orders of magnitudes lower than the minimum 
X-ray luminosity of the FGs in our sample. Thus, we 
conclude that AGN and starburst contamination does 
not contribute significantly to the X-ray measurements, 
except for AGN that might be detected only in X-rays. 
However, as noted bv lAnderson et ahl (j2007D . these ac- 
count for only a minor fraction of all AGN. 

2.2.3. How extended is the X-ray emission? 

An important issue to consider is whether the observed 
X-ray luminosity comes from an extended source. Al- 
though A GN were excised from o ur sample following the 
recipes in iSorrentino et all (l2006i), some AGN might only 
be detected in X-rav (|Tozzi et al.l 120061 : lAnderson et al.l 
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Fig. 7. — Distribution of FWHM values for 1232 point sources in 
the RASS Bright Source Catalog. The vertical dashed line marks 
the position of the peak of the distribution. The location of the 
peak, FWHM, and the corresponding la error are shown in the 
upper part of the figure. 



|2007| ). representing a further source of contamination in 
the FG sample. To objectively establish the extended 
nature of the X-ray emission we proceed as follows: 

1 - We select 1 232 point sources fr om the RASS Bright 
Source Catalog (|Voges et al.l |1999() , requiring that the 
sources have a high probability of being a real detection 
(i.e. source detection likelihood > 100) and a source ex- 
tension of zero. This data set was used to measure the 
mean FWHM of the RASS PSF over the whole sky. We 
used 2DPH0T (La Barbera et al 2008) to fit a Moffat 
distribution to every source, after a gaussian smoothing 
(with a CT of one pixel) of the original count rate image. 
The smoothing was done similarly for the FGs (see be- 
low) . Fig. [7] shows the distribution of the FWHM for 
the 1232 point sources. The location of the peak of the 
distribution (and its error) was ob tained from the bi- 
weight estimator ()Beers et al.lll990D and is found to be 
2.66(±0.02) arcmin. The narrowness of the distribution 
confirms the expectations that the RASS PSF has little 
off-axis variation resulting from the scanning strategy. 
We also verified that fitting the unsmoothed images gives 
a mean FWHM of 2.15 arc min, in good agreement with 
the value of 2.1 arcmin from lde Grandi et al.l ()1997| ). 

2 - For the seventy-four FGs, we smoothed the orig- 
inal RASS images with a 1-pixel Gaussian (as for the 
point sources). Source detection was performed using S- 
Extractor with a detection threshold of 2a over an area 
of 5 pixels, using a tophat detection filter, which is opti- 
mal for faint source detection. The gaussian smoothing 
enhances the contrast over the background. For each 
FG, we identified the closest X-ray source to the SDSS 
position. 

3 - Using the smoothed X-ray images, we fit a Moffat 
distribution to the FG candidates. Here, the error in the 
FWHM was estimated from one-hundred Monte Carlo 
simulations. For each simulation, all the detected sources 
were masked from the original RASS image and the re- 
maining pixels were bootstrapped to create a random 
background image, which was then gaussian smoothed. 
The Moffat model was added to this image, and the 
fitting was repeated. We subtracted in quadrature the 
mean FWHM of the RASS smoothed PSF from the 
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Fig. 8. — Extension parameter versus FWHM for the 43 FG 
candidates (filled circles) with an associated X-ray source. Empty 
circles represent the sample of twenty-two field galaxies defined in 
Sec. U] The error bars show la standard uncertainties. Sources are 
defined as extended if their error bar does not cross the horizontal 
dashed line corresponding to a zero extension. 

FWHM of each source. This quantity was defined as 
the source extension, with its error computed from the 
uncertainty on the FWHM. 

4 - We matched X-ray and optical sources by consider- 
ing only cases where the distance between the optical and 
X-ray FG positions were smaller than the FWHM of the 
X-Ray source. Only fourty-three (out of seventy-four) 
FG candidates have X-ray counterparts meeting this re- 
quirement. 

5 - We classified the FG X-ray sources as extended if 
the extension parameter was greater than at the 2(7 
level (see Fig. [S]). This selection removes eight point- 
like FG candidates, leading to a sample of thirty- five 
extended FGs. In the following section, we discuss the 
effect of SDSS spectroscopic incompleteness on this sam- 
ple. 

2.2.4. Spectroscopic incompleteness and the final sample 

The SDSS spectroscopy is incomplete in high-density 
regions due to limitations on fiber placement. We ad- 
dress this spectroscopic incompleteness issue as follows: 
for each of the thirty-five FG candidates, we select all 
galaxies with SDSS photometry fulfilling the -Dmax cri- 
terion and with magnitude between m and m -I- 1.75, 
where m is the FG r-band Petrosian magnitude. These 
are the galaxies that could possibly invalidate the fos- 
silness definition. Then, we queried the SDSS DR6 and 
NED for spectroscopic redshifts of these possible con- 
taminants. We found that for nine (out of thirty- five) 
FG candidates, there is at least one galaxy invalidat- 
ing the A(czmax) criterium, hence, disqualifying the sys- 
tem. Using photometric redshifts, we also rejected one 
system with three possible contaminants having concor- 
dant redshifts within 2a. Applying the same analysis to 
the eight point-like FG candidates (see previous section) 
leaves four such objects having no possible gap contam- 
inants. 

After this procedure, our final sample consists of 
twenty-five systems^^, which we use to characterize the 
properties of FGs. This is the largest sample of FGs 

^ ^ We note that almost all FGs in our sample have a total X-ray 
luminosity above the (bolometric) limit of 4.4 X 10*'^ ergs/s adopted 



available in the low redshift regime {z < 0.1). In Fig. [9l 
we show the X-ray contour plots for these FGs with the 
optical SDSS position overlaid. As we can see, usually, 
the centroid of the X-ray component coincides with the 
optical counterpart. 

2.3. Archive data and external comparison 

We searched the XMM, Chandra, ASCA and Einstein 
archives for any publicly available data, taken with any 
of the instruments, within a 1.5h^^ Mpc radius of the fi- 
nal samples of FGs and field galaxies (defined in Sec. [3]). 
We found data in the vicinity of only a single FG in 
the Chandra archive. None of the other archives con- 
tained any observations including objects in our final 
sample. We conclude that there are no useful archival 
data deeper than the RASS to systematically examine 
the X-ray properties of our FG and field galaxies, or ex- 
clude them from the sample as X-ray AGN. 

We also looked for FGs defined in the literature that 
are in the same redshift interval as our sample. Not sur- 
prisingly, we found only three s ystems well studied insey- 
eral contributions (e.g. JO03. iKhosroshahi et al.|[2007f) . 
although non e of them have data in the SDSS. Moreover, 
ISantos et al.l ([2007), defining a sample of thirty-four FGs 
from SDSS-DR4, find only four systems in the same red- 
shift regime as we consider here. One of them is in our 
sample of FG candidates, but with an X-ray component 
indistinguishable from the RASS PSF. The other three 
systems are not in our FG sample because the absolute 
magnitude limit of our SDSS catalog prevents finding a 
second-rank galaxy for the main targets. 

2.4. The Space Density of FGs 

The space density of the twenty-five FGs was estimated 
using the 1/Vmax statistics suggested by Avni & Bahcall 
(1980). The sky area is that of the DR4, 4783 square 
degrees. Applying various X-ray luminosity limits we 
find different integrated space densities, which we com- 
pared to those of previous studies. Vikhlinin et al. (1999, 
hereafter V99) studying an overluminous elliptical galaxy 
(OLEG) sample found a space density of 1.5 xlO^^/ifs 
Mpc~3 for Lx > 0.44 x lO^^/ifg^ erg/s while we find 3.4 
xlO^^ft.75 Mpc^"^. For the same limiting value of Lx, 
Romer et al. (2000) find 6.7 x 10~^fe^^ Mpc~^ JO03 finds 
6.75 xlO^'^ft.fg Mpc-^ and iDariush et al.l (|2007f) pre- 
dicts 5.4 xlO^^/iyg Mpc~'^ based on the Millennium sim- 
ulation. Considering the errors in these measurements, 
as presented in Table 1 of JO03, we can conclude that 
all these values are consistent. When a brighter limiting 
Lx is considered, Lx > 0.89 x lO^^/j-^ erg/s, V99 finds 
8.1 xl0-*^/i?5 Mpc-3 and we find 2.83 xlO'^h^ Mpc"^ 
still consistent. We emphasize that the samples analyzed 
by ROO, V99, and JO03 have only 3, 4, and 5 systems, 
respectively, while we have twenty-five new FGs, yet the 
space densities are still very consistent. 

3. COMPARING THE PROPERTIES OF FGS AND FIELD 
GALAXIES 

To establish a benchmark for the properties of FGs, 
we generated a control sample of galaxies with the same 

by Jones et al.(2003). Only two FGs (FG 11 and FG 21, see Tab. 2) 
have an X-ray luminosity of ~ 2 X 10*'^ ergs/s, corresponding to a 
bolometric value of ~ 3.5 X 10*^ ergs/s. 
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Fig. 9. — X-ray emission contours from the RASS for the twenty- 
five FGs, with the position of the optical SDSS galaxy image 
marked by a black cross. The contours have been obtained by ap- 
plying a 2 pixel gaussian smoothing to the original RASS images 
and are equally spaced at intervals of la starting at 0.5a above the 
background. The reddest contours mark the centroid of the X-ray 
emission. Each panel is three times as wide as the FWHM of the 
X-ray source. 



luminosity distribution and the same optical and X-ray 
properties as the FG sample. The control sample repre- 
sents an idealization of a sample of systems which were 
probably formed by a single collapse of a protogalaxy 
and did not experienced any major mergers in last 3-4 
Gyrs. From now on, we refer to the control sample as 
'field' galaxies. This was created by first excluding the 
brightest galaxies of the FG optical candidates and their 
companions (see Sec. 12. 1[) from the whole galaxy cata- 
log, yielding a reduced catalog with 90259 targets. We 
binned the sample of 578 FG optical candidates and the 
reduced galaxy catalog by r-band absolute magnitude, 
and for each bin computed the ratio between the number 
A^i of FG optical candidates in that bin to the number 
A'2 of galaxies in the reduced catalog. Galaxy absolute 
magnitudes were computed as described in App. |B1 To 
obtain a sample of 'field' galaxies with the same LF as 
that of the seed galaxies of the 578 optical FG candidates, 
we randomly extracted galaxies from the reduced cata- 
log in each magnitude bin according to the corresponding 
value of the ratio / = Ni/N2. Because we have excluded 
the brightest galaxies and their companion galaxies in 
the optical FG candidates from the entire galaxy cata- 
log, it is possible for N2 to be lower than A^i, yielding 
an ill-defined value of /. This is shown in Fig. [101 where 
we compare the LF of the brightest galaxies in the op- 
tical FG candidates with that of galaxies in the reduced 
catalog. At magnitudes brighter than Mr ^ —23, the 
number of optical FG candidates in each magnitude bin 
is slightly larger than the number of galaxies in the re- 
duced catalog. To ensure that the value of / is always 
less than one, we enlarged the size of the brightest mag- 
nitude bins such that the condition N2 < Ni was always 
fulfilled. This goal was achieved by considering the same 
magnitudes bins as in Fig. [10] for Mr > — 22.8, and just 
one bin at brighter magnitudes, from Mr = —23.6 to 
Mr ^ —22.8. The list of field galaxies was extracted ac- 
cording to this binning procedure, and, by construction. 



it has the same sample size (N = 578) as the optical FG 
candidates. Fig. [TOl which also plots the LF of this field 
sample, shows that its LF actually matches the LF of the 
seed galaxies in the optical FG candidates. As expected 
from the fact that we use a single, larger magnitude bin 
at ^ —23, there is only a small difference between 
the two LFs in this magnitude range. However, this dif- 
ference is only marginally significant, provided that un- 
certainties on number counts are taken into account. 

To perform a detailed comparison between FGs and 
field galaxies, we extracted a subsample of field galax- 
ies, following the same procedure adopted to select the 
sample of FGs from the initial list of 578 optical FG 
candidates (see Sec. 12.21) . We removed spirals and AGN 
from the catalog of field galaxies, and analyzed the X-ray 
images of the remaining systems. Sixty-six (out of 578) 
field galaxies have significant X-ray detections. From 
these, following the prescriptions outlined in the previous 
section, we kept only the twenty-two galaxies associated 
with an X-ray source and with a projected distance from 
rich Abell clusters larger than 1.5Mpc. Applying the X- 
ray extent requirement leaves a final sample of seventeen 
field galaxies, shown in Fig. [8] A similar analysis to that 
performed in Sec. 2.2.2 for FGs shows that contamina- 
tion from elusive AGN and starbursts is also negligible 
for the sample of field galaxies (see Fig. 6). Although se- 
lected in the same way from an optical catalog with the 
same number of galaxies (578) and the same galaxy LF 
as the FGs, the number of galaxies in the field sample 
with significant X-ray detections {N — 66) is 65% less 
than the corresponding FG sample {N ~ 102). The final 
FG and field samples have a similar ratio of 68%. 

In the following sections, we compare the distributions 
of several properties (luminosity, structural parameters, 
internal color gradients, stellar populations) of FG and 
field galaxies. For each quantity, the comparison is per- 
formed (i) by a Kolmogorov-Smirnov (hereafter KS) test 
and (ii) by deriving the mean and width of each distri- 
bution. The mean and width values obtained for each 
quantity and the results of the KS tests are reported in 
Tab.[T] In order to reduce the effects of outliers in each 
distribution, the corresponding mean, and width values 
were computed using a 2.5f7 clipping method. We also 
verified that computing the location and width values 
with ROSTAT (Beers et al. 1990) does not affect our 
results. All quantities used in this analysis are reported 
in Tab. 2 for each FG and field galaxy (see App. A). Al- 
though not used in the present analysis, the table also 
lists the properties of the four FGs with point-like X-ray 
emission that survived the final section step described in 
Sec. 12.2.41 The distributions of photometric and X-ray 
properties of FGs and field galaxies (hereafter also re- 
ferred as the field sample, FS) are presented in Fig. [11] 
and described in the following sections. 

4. LUMINOSITIES OF FIRST AND SECOND RANK 
GALAXIES 

In this section, we analyze the properties of FGs by 
examining the magnitude distributions of their seed and 
second rank galaxies as well as their X-ray fluxes. 

Fig. Illb compares the optical luminosity distributions 
of FG first rank galaxies with those of the corresponding 
field sample. The mean and width values of each distri- 
bution are reported in Tab. [T] Both the figure and the 
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TABLE 1 

Mean and width values of the distributions of several quantities for FG and FS 
GALAXIES. The p-values from the KS test comparing the distributions of FG and FS 
properties are reported in Col. (6). Values smaller than 0.1 indicate that the two 
samples likely do not derive from the same parent distribution. 





FG 




FS 




KS 




mean 


width 


mean 


width 




(1) 


(2) 


(3) 1 


(4) 


(5) 1 


(6) 






-22.64 ±0.09 


0.37 ±0.07 


-22.45 ±0.09 


0.37 ±0.05 


0.36 


5n 


2.5 ±0.4 


1.56 ± 0.30 


2.5 ±0.5 


2.0 ±0.4 


0.99 


log Rc (kpc) 


1.13 ±0.08 


0.29 ± 0.06 


0.97 ±0.08 


0.35 ±0.10 


0.23 


n 


5.4 ±0.6 


2.4 ±0.6 


5.5 ±0.6 


2.5 ±0.4 


0.83 


04 * 100 


0.11 ±0.26 


0.6 ± 0.3 


0.3 ±0.2 


0.75 ±0.11 


1.00 


V(g - r) 


-0.072 ±0.026 


0.09 ± 0.06 


-0.17 ±0.08 


0.25 ±0.08 


0.74 


logLxC-lO^^hfs^ergs/s) 


-0.90 ±0.09 


0.41 ± 0.13 


-0.84 ±0.15 


0.52 ±0.11 


0.44 




0.0 ±0.3 


0.8 ± 0.4 


-0.5 ±0.3 


1.1±0.2 


0.29 


log (TO {km/s) 


2.43 ±0.01 


0.07 ±0.01 


2.42 ±0.04 


0.09 ±0.05 


0.40 


log Age (Gyr) 


0.72 ±0.03 


0.14 ±0.03 


0.77 ±0.04 


0.16 ±0.04 


0.10 


[Z/H] 


0.44 ±0.03 


0.17 ±0.02 


0.34 ±0.06 


0.24 ±0.05 


0.29 


[a/Fe\ 


0.27 ±0.02 


0.10 ±0.01 


0.28 ±0.03 


0.12 ±0.02 


0.93 
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Fig. 10. — LFs of galaxies in the reduced catalog (empty circles), 
brightest galaxies in FG candidates (filled circles) and of galaxies 
in the field sample (diamonds). The symbols used for each sample 
are shown in the upper left part of the figure. The error bars 
mark la Poissonian uncertainties on number counts of FG brightest 
galaxies. 



table show that the r-band LFs of FG and field galax- 
ies are indistinguishable. This result is also confirmed 
by the KS test (see Tab. [T]), whose p-value is > 0.3. 
Fig. [TTb shows that not only the optical but also the X- 
ray luminosities of FG galaxies are consistent with those 
of field galaxies. This is also clear from the mean, width, 
and KS test results reported in Tab[TJ The fact that FG 
and field galaxies have the same X-ray luminosities is 
of particular interest, since it implies that both samples 
consist of galaxies embedded within halos of the same 
mass. Hence, any difference between FG and field sam- 
ples is not related to the selection of these samples, but 
is more deeply related either to the initial conditions of 
the galaxy LF within the halos, or to the different evo- 
lutionary history of galaxies within the halo. 

To analyze the luminosity distribution of second rank 
galaxies, we have to consider that the SDSS-DR4 galaxy 




-2 2 
□ ^(.100) 

Fig. 11. — Distributions of FG and FS properties as follows: (a) 
absolute magnitude, (b) X-ray luminosity, (c) density excess, (d) 
normalized distance to the red sequence, (e) r-band effective radius, 
(f) r-band Sersic index, (g) internal color gradient, (h) boxiness 
parameter. FG and field galaxies are plotted as filled and empty 
circles, respectively (see upper— left panel). Error bars mark the Icr 
standard confidence intervals of FG number counts. To make the 
plot more readable, uncertainties on field galaxy number counts 
are not shown. The distribution of field galaxies is normalized to 
the total number of corresponding FG galaxies. 



catalog is magnitude limited and that FG optical candi- 
dates are selected as systems where the magnitude differ- 
ence between the two brightest galaxies is greater than 
AAfniin- This implies that a larger fraction of FGs can 
be selected at the brightest magnitudes of the galaxy 
LF. In other words, the seed galaxies of FGs do not form 
a magnitude complete sample. Such incompleteness is 
illustrated in Fig. [121 where we plot the magnitude dif- 
ference AM between the second and first rank galaxies 
of FGs as a function of the r-band absolute magnitude, 
"■^Mr, of the first rank galaxy. Absolute magnitudes 
are obtained as described in App. [B] by k-correcting 
the Petrosian r-band magnitudes of the DR4 dataset 
to redshift z ~ 0.1. Since the SDSS galaxy catalog is 
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limited to ~ —20, and since, by definition, we have 
AM >AMmin= f .75mag, the seed galaxies of FGs have 
to fulfill the equation 

O-^Mr +°-^ kr < -20 - AM 21.75mag. (1) 

The term ^'^k^ denotes the k-correction computed as de- 
scribed in App.IbI and is introduced in the previous equa- 
tion to account for the fact that k-corrections were not 
applied in the computation of the r-band absolute mag- 
nitudes used to select the SDSS-DR4 catalog. Since the 
typical value of k^ is — O.lSmag, Eq. [1] implies the con- 
straint °-^Mi- < —21.62, and in fact, as shown in Fig. [T^ 
we do not find any first rank galaxy fainter than this mag- 
nitude limit. Moreover, for a given value of '^'^M^, we can 
only select FGs with °-^M, kr + AM < -20, which 

implies AM < -°-'^M, - 20 -O ^ k,. "-^M, - 19.87. 

This last constraint is shown by the line in Fig. [121 and 
results in a larger range of allowed values of AM at 
bright galaxy magnitudes. In other words, the selection 
procedure affects the magnitude distribution of FGs, bi- 
asing the sample towards a larger allowed fraction of sys- 
tems at brighter magnitudes. Despite this bias. Fig. [T^] 
shows that the brightest first rank galaxies do not show 
a larger range of AM values, although they are allowed. 
At bright magnitudes, we see that seed galaxies of FGs 
are characterized by a larger magnitude gap. However, 
one has to note that the SDSS spectroscopic incomplete- 
ness due to fiber collisions might cause an overestimate 
of AM, and that the SDSS spectroscopic survey did not 
target many very bright galaxies (r < 15.5), which might 
reduce the number of detected FG systems at the bright 
end. 

To better understand the distribution of AM with 
respect to the magnitude of the seed galaxies, we ran 
Monte-Carlo simulations. For each seed galaxy of the 
102 FG candidates, we assigned a fictitious second rank 
galaxy from the SDSS catalog by randomly extracting 
a galaxy with AM > AMmin. We repeated this proce- 
dure several times {N = 10000), and estimated the mean 
value of AM obtained for the randomized second rank 
galaxies. Fig. 1121 shows that this mean AAI value actu- 
ally decreases as '^■^M^ increases. At magnitudes fainter 
than '^■^Mr ^ —23, the simulated AM value is fully con- 
sistent with the AM distribution of FGs, implying that 
the AM vs. °-^Mi- correlation can be simply explained 
by the shape of the LF of second rank galaxies. However, 
at bright magnitudes {^-^M^ < —23), all the seed galax- 
ies of FGs have AM values larger than the simulations, 
which might reflect some intrinsic physical property of 
FG systems (see Sec. (HI)- 

5. FAINTEST GALAXIES: DENSITY EXCESS 

We characterize the population of the faintest galaxies 
around each FG by measuring the density excess of these 
galaxies with respect to the field population. Since galax- 
ies around FGs can be fainter than the magnitude limit 
of the SDSS spectroscopic catalog, we measure the den- 
sity excess by taking advantage of the SDSS-DR4 pho- 
tometric d ata, whose completeness limit^-^ is ~ 22.2mag 
in r-band (jSthoughton et aDl2002l ). The density excess, 
(5n, is defined as follows: 



This is defined as 95% detectabiUty for point-like sources. 




Fig. 12. — Differences between r-band magnitudes of second 
and first rank galaxies of FGs, AM, are plotted as a function of the 
absolute magnitude, ''■^Afr, of the first rank galaxies. Grey circles 
denote the sample of 102 FG candidates with X-ray significant 
detection (see Sec. 12.21 1. while black symbols correspond to the 
final sample of twenty-five FGs. The black line marks the range of 
allowed values of AM as a function of "'^Mr, while the grey curve 
shows the mean AM value as a function of ^'^M^, obtained by 
randomly extracting second rank galaxies from the SDSS catalog 
(see the text). 



where p-m is the mean galaxy density within a inner cir- 
cle of radius rjn centered on the FG seed galaxy and 
Pout is the mean galaxy density in a concentric annu- 
lus with inner and outer radii of ri and r2, respectively. 
The quantities din and (Tout denote the Icr uncertainties 
on Pin and pout, which are estimated by accounting for 
Poisson noise in the galaxy counts. For each FG, we 
measure (5n by setting ri„ = 0.25i?Abcib — 1-RAbcii 
and r2 — 2RAbeU, respectively, where RAbeii is the ap- 
parent size corresponding to one Abell radius (1.5Mpc) 
at the redshift of the FG seed galaxy. We have verified 
that our results do not change when these radii are var- 
ied by as much as 20%. To probe a fixed range of the 
galaxy LF, we measure pin and Pout using only galaxies 
in the magnitude range of m2 to m2 -t- 2, where TO2 is the 
magnitude of the second rank galaxy in each FG. 

Fig. [Tib compares the 5^ distributions of FGs and field 
galaxies. Since we do not have second rank galaxies for 
the field sample, we obtained the corresponding density 
excesses as follows. For each field galaxy with apparent 
and absolute magnitudes and "'^M^, we randomly 
assign a magnitude shift AM following the distribution 
of magnitude differences between first and second rank 
galaxies of the 102 FG candidates (see Fig. [T2)) . Then, 
we set m2 = fn^ -\- AM, and measure the density ex- 
cess around the field galaxy following the procedure out- 
lined above. In other words, we assign a fictitious second 
rank galaxy to each field galaxy according to the magni- 
tude distribution of first and second rank galaxies in the 
sample of 102 FG candidates. Fig. [TTb shows that most 
galaxies in each sample tend to have positive 5^ values. 
In fact, as shown in Tab.[l] the mean value of the density 
excess is about 2.5 for both samples, and, considering the 
corresponding uncertainty, is always significantly larger 
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than zero. This imphes that there is an excess of faint 
galaxies in the surrounding regions of both FG and field 
galaxies. Both Fig. [Tib and the mean and width val- 
ues reported in Tab. [1] show that the distributions of 5n 
values of FG and field galaxies are fully consistent. 

6. FGS AND THE RED SEQUENCE 

The distribution of FG seed galaxies in color- 
magnitude space can provide clues to their formation 
history. We characterize the g — r vs. r color-magnitude 
relation of early-type galaxies as described in App. [b1 
and then compare the distribution of FG and field galax- 
ies to that relation. It is implicit in this simple analysis 
that the red sequence mainl y originates from a mass- 
metallicity relation (see e.g. iKodama et anil998D . and 
that galaxies lying on the sequence represent, as a first 
approximation, passively evolving systems that have not 
experienced recent episodes of star formation. By mea- 
suring how distant the colors of FGs are from the color- 
magnitude relation, we have an indication of how star 
formation is evolving in these systems. 

We compute the distance 5cm of the colors of FG and 
field galaxies to the red sequence as 



(3) 



where a and b are the offset and slope of the color- 
magnitude relation, and C o-rC "'"'" Afr) is the scatter about 
the red sequence (see Eq. lBll of Appendix [B|) . Fig. fTTH 
compares the Scm distribution of FG and field samples, 
with the corresponding mean and width values tabulated 
in Tab.[T] The results of the KS test, comparing the Scm 
distributions, is also reported in the same table. We see 
that the 6cm distributions are always peaked near zero, 
as one would expect since FGs and field galaxies are mor- 
phologically selected as early-type systems. Fig.fTTH sug- 
gests that field galaxies have a broader distribution than 
FGs. However, the corresponding mean values are fully 
consistent with zero and their widths are consistent with 
unity, implying no significant difference between the Scm 
distributions. This result is also confirmed by the p- value 
from the KS test. 

7. STRUCTURAL PROPERTIES AND COLOR GRADIENTS 

We analyze the surface brightness distribution of FG 
seed galaxies by studying (i) their structural parame- 
ters and internal color gradients, and (ii) the departures 
of galaxy isophotes from purely elliptical shapes. Inter- 
nal color gradients are estimated using the structural pa- 
rameters of galaxies, specifically the effective radius, r^, 
and the Sersic index (shape parameter), n, while galaxy 
isophotes are characterized using the 04 parameter that 
describes boxy /disky departures of the isop hotes from el- 
liptical shapes ()Bender fc Mollenhoff ]|1987D . The deriva- 
tion of these quantities is described in App. ICl 

Figs.fTTb andfTTF compare the distributions of effective 
radii and shape parameters for FG and field galaxies. 
The mean and width values of FGs and the field sam- 
ple (see Tab. [1]) are fully consistent, with the p-values 
of the KS tests greater than ~ 0.2. The lack of galaxies 
with Sersic index n > 2 is a natural consequence of our 
sample consisting of morphologically selected early-type 
galaxies. Although the Sersic index does not provide a 
precise measurement of the bulge to disk ratio of a galaxy, 



it can effectively separate disk and bulge dominated sys- 
tems, with n ~ 2 as the dividin g value between the two 
classes (e.g. iBlanton et al]|2005[ ). 

Fig. [TTk compares the distributions of color gradients, 
V{g ~ r), for FG and field galaxies. We find that fossils 
and field galaxies have a consistent mean color gradient of 
V{g~r) ~ —0.07, but the scatter of field galaxies around 
the mean value is somewhat larger than that of FGs (see 
Tab.[T]). The difference between the width values of the 
two samples is 0.16 ± 0.1, which is (marginally) greater 
than zero by 1.6cr. The KS test (Tab. [l]) confirms that 
there is no significant difference between the two sam- 
ples. Our color grad ient estimates can be dir ectly com- 
pared with those of iLa Barbera et all (|2005D . who de- 
rived the mean g — r color gradient in early-type galaxies 
in clusters of different richnesses at different redshifts. 
For galaxies in poor groups at a redshift of z ~ 0.08, 
close to the median redshift of our FG and FS samples, 
they found a mean color gradient of -0.0641^ °^ (see 
their Tab. 1), agreeing with the mean V(g — r) values 
reported in Tab. [TJ The mean val ue of V(q — r) ha s also 
been measured from SDSS data by IWu et all (|2005D , who 
analyzed a sample of 36 early-type galaxies at z ~ 0.02. 
They find a mean V{g — r) value of — 0.05±0.01 (see their 
Tab. 3), which differs by less than 2a from the results in 
Tab.IH 

In Fig. Illh . we compare the distributions of 04 val- 
ues for FG and field galaxies. Looking at the mean and 
width values reported in Tab. [T] as well as at the KS test 
results reported in Tab. [1] the distributions for FG, and 
field galaxies are again fully consistent. The distribution 
of 04 values carries interesting clues about the formation 
of FG galaxies. Analyzing the isophotal shapes of seven 
FG galaxies drawn from the sample of Jones et al. (2003), 
Khosroshahi et al. (2006, hereafter KPJ06) found that all 
FG galaxies tend to have disky isophotes. In more de- 
tail, examining Fig. 3, three (of seven) FGs in the KPJ06 
sample have 04 values very close to zero, with the re- 
maining four galaxies having significantly larger values. 
In our sample, we find that the distribution of FGs is 
peaked around 04 = 0, with equal numbers of disky and 
boxy systems. The disagreement we find with KPJ06 
may reflect the small sample analyzed in that study. Our 
comparison of FG and field galaxies shows no difference 
in isophotal deviations from elliptical shapes. In other 
words, FG galaxies do not have peculiar isophotes when 
compared to 'normal' field galaxies. 

8. STELLAR POPULATIONS 

To analyze the stellar populations of FGs and field 
galaxies, we use spectra from the SDSS-DR4 with a res- 
olution of 4.42 A (FWHM). For each galaxy, we search 
for the Single Stellar Population (SSP) model that best 
matches the galaxy spectrum. The matching is per- 
formed using measurements of spectral indices, resulting 
in estimates of age, metallicity, and elemental abundance 
ratios [a/Fe]. The procedure consists of (i) estimating 
the spectral indices, and (ii) extracting the stellar popu- 
lation parameters, as detailed in App. |DJ 

Velocity dispersions and stellar population parameters 
for FGs and field galaxies are reported in Tab. 2. The 
distributions of cr, age, [Z/H], and [a/Fe] are shown in 
Fig. 1131 with the corresponding mean and width values, 
and the KS test being reported in Tab. [T] We see that 
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Fig. 13. — Same as Fig. Illl but comparing velocity dispersions 
and stellar population parameters for FG and field galaxies. 
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Fig. 14. — Correlations of stellar population parameters with 
velocity dispersion for FGs, field galaxies, and galaxies in Hickson 
Compact Groups, represented with different symbols as shown from 
the inset in the middle panel. In the bottom figure, as an aid to 
the discussion, the best fit to the sample of sixty-six field galaxies 
(see Sec. |3J is plotted (dashed line). Values of stellar population 
parameters from the de la Rosa et al. (2007) have been transformed 
to the a-models system used here (see App.[D}. 



the distributions of FGs and field samples are fully con- 
sistent, with the p- value of the KS test being larger than 
~ 0.1. Fig. [13] shows how stellar population parame- 
ters vary as a function of velocity dispersion for both 
FGs and field galaxies. We see that both samples follow 
the same correlation between [a/Fe] and velocity disper- 
sion, as found from previous studies for field galaxies (e.g. 
de la Rosa et al. 2007). Moreover, both FGs and field 
galaxies occupy the same region in the diagrams of age 
and metallicity versus velocity dispersion, emphasizing 
the similarity of both samples. In Fig. [TH we also insert 
the stellar population parameters for the twenty elliptical 



galaxies in Hickson compact groups (HCGs) analyzed by 
de la Rosa et al. (2007). The stellar population parame- 
ters of these galaxies were remeasured following the same 
procedure as for FGs and field galaxies. Note that ellip- 
ticals in HCGs have higher ages and lower metallicities 
than both FGs and field galaxies, showing no correlation 
between [a/Fe] and velocity dispersion. 

9. DISCUSSION 

In the current hierarchical model of galaxy formation, 
the basic idea is that massive galaxies form later as the 
result of the merging history of their dark matter halos, 
but their stars form earlier (downsizing). In this sense 
special attention is being focused on massive (luminous) 
galactic systems, whether they reside in low or high den- 
sity locations. We may find fossils at different cosmic 
epochs, and they provide important clues that help us 
distinguish between nature and nurture. Our main re- 
sults from this paper are: 

We found 578 candidates following our optical photo- 
metric definition of FGs in the redshift range 0.05 < z 
< 0.095. After discarding incorrect morphological clas- 
sifications and AGN we are left with 102 FGs for which 
there is a significant X-ray detection in ROSAT (maxi- 
mum total flux of 10~"'^^erg cm~^s~^). Only twenty-five 
FGs were selected with likely extended X-ray emission 
and not associated with rich clusters. These constitute 
our final sample of FGs. Applying the V/Vmax test we 
find a spatial density of 2.83 xlO^^hys^ Mpc"^, com- 
parable to three other independent observational studies 
(V99, ROOO, and JO03) and to the recent analysis of the 
Millennium Simulation by Dariush et al. (2007). 

As shown in Figure 12, the brighter FGs always have 
AM greater than ^ 3.0 mag while fainter FGs are limited 
by the completeness of the SDSS catalog. The lack of 
systems in the lower left part of the figure, at magnitudes 
brighter than °-^Mr ^ —23, might be caused by selection 
effects as well as by the shape of the galaxy LF. Thus, 
it is not clear from the present analysis if the absence 
of systems in the lower left part of Fig. 12 is related to 
some physical characteristics of the sample. It would be 
reasonable to expect that as mergers progress the LF of 
such systems will have a larger gap between the first and 
the second ranked galaxies. 

The density excess around FGs is fully consistent with 
what we measure for field galaxies, implying that the 
environments around both systems are statistically iden- 
tical when considering the faint galaxy population. This 
result seems to indicate that FGs and field galaxies orig- 
inate from similarly overdense regions of the Large Scale 
Structure. 

By studying the distribution of color offsets from the 
red sequence, we find that FGs occupy a similar region 
to field ellipticals. We also find that the mean color gra- 
dients in both samples (FG and FS) are statistically the 
same. However, the scatter around the mean for FGs is 
significantly (1.6 a) smaller than that for field ellipticals, 
which may indicate a more regular process in the buildup 
of FGs, such as mergers of L* galaxies. The distribution 
of a4, which measures isophotal shape deviations from a 
pure ellipses, indicates that FGs and field ellipticals are 
similar. 

Studying the stellar populations in the twenty-five FGs 
and seventeen field galaxies we find that there is no dif- 
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ference in age, metallicity or a-enhancement, indicating 
that the star formation history of fossil groups seems to 
be analogous to that of field ellipticals. We further ex- 
amined elliptical galaxies in compact groups studied by 
de la Rosa et al. (2007). As already established in previ- 
ous works, elliptical galaxies in CGs are older and more 
metal-poor than field ellipticals (Proctor et al. 2004; de 
la Rosa et al. 2007) and fossil groups. Particularly, 
when we look at the relation between [Z/H] and logcr, the 
mass-metallicity relation is evident (see Lee et al. 2006; 
Kobayashi 2005). More massive galaxies (in both FGs 
and field galaxies) retain their metals more efficiently 
than the low mass ones (i.e., ellipticals in CGs). The a- 
enhancement diagram shows a clear trend between a/i^e 
and velocity dispersion (mass). This can be interpreted 
as a manifestation of downsizing - the less massive galax- 
ies have a more extended star formation history. Based 
on these results we see that FGs are similar to field ellip- 
ticals but cannot be formed by dry mergers of ellipticals 
in a CG. However, a wet merger with a gas rich disk 
system may also explain the observed relations. Here, 
the luminosity- weighted parameters (Age, [Z/H], [a/Fe]) 
are more sensitive to the recent star-formed-population, 
generating lower ages, higher [Z/H] and higher [Mg/Fc]. 
This assumes that the merger-starburst is soon followed 
by SN/AGN feedback which depletes the gas content. 

McCarthy et al. (2004) have found a significant pop- 
ulation of red galaxies at high redshift (1.3 < z < 2.2), 
and claim that they cannot be descendants of the z ^ 3 
Lyman-break galaxies. Based on their spectroscopic data 
they also conclude that most of the present day mas- 
sive galaxies had an early (zf — 2.4) and rapid forma- 



tion. If we restrict ourselves to the twenty-five FGs, their 
last main star formation episode must have occurred at 
z > 0.3, much more recently than in the high-z galaxies. 
It appears that a more recent burst occurred in the FGs 
which would explain their lower ages. Therefore, we may 
be seeing massive ellipticals that accreted a small galaxy 
at z > 0.3. 

Our results indicate that FGs are not significantly dif- 
ferent from bright ellipticals found in low-density regions 
of the universe. As claimed by Mulchaey & Zabludoff 
(1999) in a detailed study of NGC 1132, FGs could sim- 
ply be "failed" groups - a few relatively bright galax- 
ies merged forming the dominant system we see today, 
but without enough surrounding matter to form addi- 
tional nearby bright galaxies, resulting in an atypical LF. 
This is also in agreement with recent results obtained by 
D'Onghia et al. (2005) based on N-body/hydrodynamical 
cosmological simulations where they find that FGs are 
formed by the infall of L ~ L» systems along filaments 
with impact parameters as small as 5 kpc. 

Our findings suggest that the objects meeting the ob- 
servational criteria expected for fossil groups at low red- 
shift might not be true fossils. The extensive study pre- 
sented here whereby we select FG candidates based on 
optical data, and then examining their X-ray counter- 
parts reveals no difference between FGs and field ellip- 
ticals. We note that fossils can be created anytime in 
cosmic history, and these systems may represent the fi- 
nal stage of mass assembly as suggested in the analysis 
of the Millennium simulation by Dariush et al. (2007), 
instead of forming a distinct class. 



APPENDIX 

LIST OF PROPERTIES OF FGS AND FIELD GALAXIES 

Tab. 2 lists all the quantities used to characterize the samples of FGs and field galaxies. FGs with extended X-ray 
contours {N = 25) have flag=l in column 21, while field galaxies {N = 17) have flag=2. The table also lists the four 
FG candidates that are not in the final FG sample because they do not fulfill the X-ray extension criterion. These 
objects have flag=0. Column 1 is a running number for the list of FGs and field galaxies. Columns 2, 3, and 4 give the 
RA and DEC (in degrees) and spectroscopic redshift from the SDSS DR4. Columns 5, 6, and 7 provide the z = 0.1 
k-corrected absolute magnitude, the magnitude gap of the FG, and the g ~r color, respectively. For field galaxies, the 
values of AM have been randomly assigned as described in Sec. ID Column 8 provides the density excess around each 
galaxy. Columns 9 and 10 report the effective radii (in arcsec) in the g and r bands, while columns 11, and 12 list 
the corresponding Sersic indices. Columns 13 and 14 provide the a4(xl00) parameters and internal color gradients 
(in mag/dex). Column 15 reports the X-ray luminosity in units of lO^^/ifg^erg/s. Column 16 provides the distance of 
each galaxy to the red sequence (Eq. [3|). Columns 17 to 20 provide the velocity dispersions (in units of km/s) and 
stellar population parameters. 

MEASURING THE RED SEQUENCE 

To quantify the properties of the FG first-rank galaxies relative to the red sequence, we use the SDSS-DR4 galaxy 
catalog to calculate the offset, sl ope and scatter of the 5 — r vs. r color-magnitude relation for all early-type galaxies. 
Following a procedure similar to iBernardi et al.l ()2003D . we select early- type galaxies on the basis of the SDSS spectro- 
scopic parameter edass, which classifies the spectral type using a principal component analysis, and the photometric 
parameter fracDeVr, which measures the fraction of the galaxy light distribution that is better fit by a de Vaucouleurs 
law. We define as early-type galaxies those systems with eclass<0 and fracDeVr >0.8, resulting in a sample of 39733 
(out of 91563) galaxies selected from the SDSS catalog used to search for fossil systems (see Sec. 12.11) . Total galaxy 
magnitudes are taken from the r-band Petrosian magnitudes output by the SDSS photometric pipeline, since they are 
model-independent and highly reliable for bright galaxies. Colors are obtained from the g and r band SDSS model 
magnitudes because they are measured in the same aperture for all filters. Both magnitudes and colors are k-corrected 
using the fccorreciv4_l_4 software package (Blanton et al. 2003), allowing k-corrections to be estimated through filters 
that are blue-shifted by a factor (1 -I- zq). This procedure has the primary advantage of minimizing uncertainties on 
k-corrections, provided that the value of zq is close to the redshift range of the observed galaxie s. Following othe r 
papers that analyze the color and magnitude distribution of galaxies in the SDSS database fe.g. iHogg et al.ll2004j) . 
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we choose zq — 0.1, corresponding to the upper redshift limit of our galaxy catalog, and we indicate the zq = 0.1 
k-corrected model magnitudes with '^■^g and "-^r, and the zq = 0.1 k-corrected Petrosian absolute magnitude with the 
notation "-^Mr. 

Fig. [15] shows the color-magnitude diagram of both early- type and all galaxies in our catalog. Although the catalog 
is almost complete down to Mr — —20 (see Sec. 12. ip . we note that the difference in k-corrcctions for galaxies with 
Mr ~ —20^^ makes the color-magnitude diagram incomplete down to that magnitude limit. Therefore, we derive the 
slope, a, and the offset, b, of the color-magnitude relation by minimizing the rms of galaxy colors around the linear 
relation, using only gal axies brighter than Mr ~ —20.2. To reduce the effect of outliers, the color rms is computed via 
bi- weight statistics (Be ers et al.lll99dt ). The uncertainties on a and b are estimated through the bootstrap method, with 
N = 1000 bootstrap iterations. For the sample of early- ty pe galaxies, we find a = 0.944±0.005 and b = — 0.0206zb0.003 , 
consistent with b = — 0.022 reported bv lHogg et al.l (|2003l ) and with other results in the literature (e.g. lCool et al.|[2006t 
iGallazzi et al.l[2006f l. The derived values of a and b change by less than 10% if we apply the same fitting procedure 
to the entire galaxy catalog (i.e. not applying the eclass and fracDeVr cuts), demonstrating the robustness of the 
regression method. To estimate the scatter around the color-magnitude relation, we bin the sample of early-type 
galaxies with respect to magnitude, and derive the color rms, cTg-r, in each bin with bi- weight statistics. Fig. 1161 plots 



as a function of ^-^Mr. We find that CTq 



relation 



(7g^r{ ' M, 



increases at the faintest magnitudes, with this trend well-fit by the 
0.045 - 0.00047 x (O ^A/^ -I- 20.2)"*, (BI) 



in the magnitud e range —23 <°'^ Mr 
studies (see e.g. lCool et al.l[2006[ ). 



< —20.2. The mean value of (Jg-r is 0.035, consistent with that found in previous 



MEASURING STRUCTURAL PARAMETERS AND INTERNAL COLOR GRADIENTS 

Structural parameters were derived using 2DPH0T (La Barbera et al. 2008, hereafter LdC08), an automated software 
environment that performs several tasks, such as catalog extraction (using S-Extractor, Bertin & Arnouts 1996), 
star/galaxy separation, and surface photometry. To estimate color gradients, structural parameters were measured in 
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Fig. 15. — Color-magnitude diagram of early-type galaxies (upper panel) and of all galaxies (lower panel) in the catalog. The grey line 
in the upper panel shows the bi- weight fit to the data, while grey circles, and error bars plot the mean and standard deviation, dg—r, of 
the color distribution in different magnitude bins. The grey solid line in the lower panel is the same as that in the upper panel, while the 
grey dotted line shows the bi-weight fit to the color-magnitude relation of all galaxies in the catalog. 



As mentioned in Sec. [S] k-corrections have not been applied when computing total magnitudes to select the SDSS-DR4 catalog 
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Fig. 16. — The scatter around the red sequence, Ug^r, plotted as a functi on o f central value of each magnitude bin in the upper panel 
of Fig. 1151 The grey curve is the polynomial fit to the points, as given in Eg. IBll 



both the g and r bands. For each galaxy in the FG and field (FS) samples we retrieved the corresponding g- and 
r-band corrected frames from the SDSS DR4 Data Archive Server and ran 2DPH0T on them. A complete description 
of the 2DPH0T package can be found in LdC08; here we provide a brief description of the steps relevant to measuring 
galaxy structural parameters. For each galaxy, a PSF model was constructed by fitting the four closest stars in the 
image with a sum of three two-dimensional Moffat functions. Isophotal distortions of star isophotes were modeled as 
described in LdC08. The parameters r^, </i>e and n were then obtained by fitting galaxy images with PSF-convolved 
Sersic models. As an example of the fitting procedure, Fig. [T7] plots the images of a subset of galaxies from the list of 
twenty-five FGs, along with the residuals after subtracting the best-fitting Sersic models. The shapes of FG and FS 
candidates are chara cterized by the isophotal sha pe parameter, 04, of each galaxy. Isophotal contours were expanded 
into a sin/cos series ()Bender fc Mollenhoff |[l987l ). and analyzed as described in LdC08, by deriving the 04 coefficient 
as the average of the 04 profile within a radial range from four times the FWHM to twice the galaxy effective radius. 

Interna l color gradients, defined as the logarithmic slope of the galaxy radial color profile, were estimated as follows 
(see also iLa Barbera et al.l [20031 ) . For each galaxy in the FG and FS samples, the corresponding g- and r- band 
best-fitting Sersic parameters were used to construct seeing-deconvolved galaxy images. A set of concentric ellipses 
were constructed for both the g- and r-band galaxy images, with all ellipses having the same position angle and axis 
ratio parameters as those derived from the r-band Sersic fitting (Sec.[C]). For each ellipse, a g — r color was calculated 
as the average ratio of the g- and r-band intensity values along that ellipse. The color profile was then computed as 
the average g — r color as a function of the equivalent radius of the corresponding ellipse. An ordinary least-squares 
fit to the color profile, adopting the color index as the independent variable, was then performed, with the slope of 
the fitted line yielding the internal color g radient for each galaxy. Following previous studies of color gradients in 
early- type galaxies fe.g. iPeletier et al]|1990r ). the fit was performed in a range of the color profile, from an inner radius 
of 0.1 to an outer radius of 1.0 times the r-band effective radius. 

SPECTRAL ANALYSIS 
Measuring spectral parameters 

The spectra of FGs and field galaxies were masked to avoid bad pixels affecting the wavelength ranges defining the 
spectral indices, brought to redshift zero and corrected for foreground dust extinction following Schlegel, Finkheiner 
& Davis 1998. Then, we derived the spectral indices H/?, [MgFe]' = (Mgb (0.72 x Fe5270 + 0.28 x Fe5335))^/^ 
Fe3 = (Fe4383-|-Fe5270-|-Fe5335)/3, and iFei= (Fe5270-hFe5335)/2, which are defined according to the Lick system 
(Worthey & Ottaviani 1997; Gonzalez 1993), except for the higher resolution of the SDSS spectra. 

The measurement of the H/3 spectral index, which is directly related to age, may be significantly affected by nebular 
emission, biasing the inferred age towards older values. To correct for this effect, we used the STARLIGHT code (Cid 
Fernandes et al. 2005) which finds the best fit to the observed spectrum by combining a set of SSP-MILES Spectral 
Energy Distributions (SEDs; see Sec. 8.2), convolved to the SDSS spectral resolution and broadened to match the 
galaxy's velocity dispersion. As a result of the fitting procedure, STARLIGHT also provides the velocity dispersion, a, 
of each galaxy. During the fitting process, regions around emission lines were masked. Fig. 18 shows how the residual 
spectrum unveils the nebular emission. Instead of correcting for the residual emission line, we measured the H/3 index 
directly from the absorption spectra obtained by STARLIGHT. 

Velocity dispersions of FGs and field galaxies were taken from the SDSS DR4. For 53 galaxies (out of 102 FG 
candidates and 66 field galaxies with significant X-ray detection), the a measurements were not available from the 
SDSS database, and in this case we adopted the STARLIGHT estimates. The spectral parameters, i.e., velocity 
dispersions and spectral indices, were corrected to a common physical aperture of Re/^/8, following J0rgensen et 
al. (1995) and J0rgensen (1997). For each galaxy, the effective radius. Re//, was defined as Reff — {b/a) x Rdev: using 
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Fig. 17. — Two-dimensional fitting of four galaxies drawn from the sample of twenty-five FGs. For each galaxy, the upper-left and right 
subpanels show the galaxy r-band image and the corresponding residual map after model subtraction, respectively. Right ascension and 
declination are reported, in degrees, in the lower right corner of the left subpanel. The spatial scale is shown in the lower-left corners. The 
lower six subpanels show the fitting residuals, normalized to the photon noise in each pixel, as a function of the radial distance, r, to the 
galaxy center within six bins of the polar angle, 9, respectively. The limits of each bin are shown in the lower right corner of each subpanel. 
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Fig. 18. — The synthetic spectrum obtained by STARLIGHT (continuous line) is the combination of SSP models which best fits the 
observed SDSS-spectrum (dotted line). The residual spectrum (observed - synthetic), can be used to measure emission indices. Instead of 
correcting for emission filling, we measure the H/3 index directly from the synthetic spectrum 
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the deVRadr (Rdev) and deVABr (b/a) parameters from the SDSS database. Since the aperture corrections rely on 
effective radii derived by fitting de Vaucouleurs models, we adopted the SDSS photometric pipeline's Reff values, 
instead of using our effective radii derived using Sersic models (see App. [Ct . 

Extraction of stellar population parameters 
We measured age, metallicity and abundance ratio [a/Fe] using the recently developed a-enhanced SSP-MILES 
models (de la Rosa et al. 2008), which provide SEDs at 2.4 A (FWHM) resolution for different a-enhancement ratios, 
combining theoretical (Coelho et al. 2005) and empirical (MILES) libraries (Sanchez-Blazquez et al. 2006). For each 
galaxy, we constructed a grid of H/? and [MgFe]' values using the solar metallicity SSP models ([a/Fe]=0) with a 
variety of age and [Z/H] combinations. The age and [Z/H] of the observed spectrum were obtained by interpolating its 
two spectral indices among those of the models. The H/3 and [MgFe]' indices are mostly sensitive to age and metallicity, 
respectively, breaking the well-known age-metallicity degeneracy. Moreover, they have the advantage of being almost 
insensitive to [a/Fe] (Korn, Maraston & Thomas 2005). As a second step, after fixing the estimated age, we construct 
a new grid of Fe3 and Mgb values with models having a variety of [Z/H] and [a/Fe] combinations. Interpolation of the 
observed indices yields the final metallicity and [a/Fe] values. To check the reliability of the a-enhancement estimates, 
we repeated the above computation, replacing the Fe3 index with <Fe>= (Fe5270 -I- Fe5335)/2. The two procedures 
are fully consistent at a confidence level of 2a for all but three of the EG and FS galaxies. For those three objects, we 
adopted the more accurate <Fe> estimate (see de la Rosa et al. 2008). 
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Table 2. Properties of FG and FS galaxies. 
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0.25 




20 


195.28142 


-3.4480 


0.083 


-22.761 


2.74 


1.007 


5.6 


10.16 


12.72 


4.9 


4.7 


0.46 


-0.086 


-1.11 


1.070 


2.425 


0.9 


0.57 


0.34 




21 


354.57353 


15.6689 


0.066 


-22.600 


2.68 


0.973 


1.7 


8.45 


14.56 


2.7 


4.0 


2.66 


0.044 


-1.77 


-0.122 


2.361 


0.8 


0.21 


0.23 




22 


129.81595 


28.8441 


0.079 


-23.281 


3.04 


0.993 


4.6 


4.51 


6.14 


2.1 


2.4 


0.20 


-0.203 


-1.10 


0.122 


2.543 


0.7 


0.78 


0.38 




23 


230.21770 


48.6607 


0.074 


-23.126 


3.15 


0.943 


4.0 


10.66 


12.67 


4.7 


4.8 


-0.72 


-0.102 


-0.40 


-3.176 


2.529 


1.0 


0.27 


0.45 




24 


228.00073 


36.6520 


0.066 


-22.827 


2.91 


0.980 


0.9 


11.04 


14.23 


3.2 


3.3 


-1.16 


0.059 


-0.49 


-0.078 


2.417 


0.7 


0.60 


0.14 




25 


234.78742 


33.6881 


0.070 


-22.555 


1.88 


1.016 


3.8 


20.37 


20.19 


9.9 


10.4 


0.09 


-0.129 


-0.84 


1 .299 


2.458 


0.7 


0.46 


0.45 




26 


172.88388 


12.6997 


0.081 


-22.306 


1.85 


0.972 


1.0 


7.00 


9.99 


2.4 


2.6 


-0.04 


-0.108 


-1.23 


0.017 


2.434 


0.4 


1.46 


0.31 






Table 2 — Continued 



ID 


R.A. 


Dec. 


z 




AM 










"g 




a4(xl00) 


V(g - r) 




ScM 


logo-Q 


log Age 


[Z/H] 


[a/Fe] 


FLA 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21; 


27 


20.09644 


-0.0790 


0.077 


-23.128 


3.01 


1.003 


1.6 


32.20 


39.10 


4.8 


5.0 


-0.42 


-0.133 


-1.00 


1.110 


2.543 


0.8 


0.40 


0.32 





28 


1 80.40768 


1.1980 


0.082 


-22.225 


2.14 


0.903 


-2.6 


10.41 


9.42 


4.7 


9.2 


2.29 


-0.848 


-0.85 


-1.773 


2.197 


0.6 


0.09 


0.03 





29 


118.18416 


45.9493 


0.052 


-22.167 


2.05 


0.963 


2.3 


6.06 


5.78 


4.5 


4.8 


-0.45 


-0.042 


-0.96 


-0.144 


2.336 


0.7 


0.23 


0.17 





30 


168.84945 


54.4441 


0.070 


-22.802 


1.99 


0.939 


2.6 


7.22 


7.38 


6.3 


5.6 


-0.89 


-0.041 


-1.42 


-1.795 


2.349 


1.0 


0.09 


0.25 


2 


31 


232.31104 


52.8640 


0.073 


-22.927 


2.17 


0.995 


3.2 


21.23 


24.43 


9.8 


8.6 


-0.34 


0.017 


-0.36 


0.602 


2.454 


0.8 


0.60 


0.41 


2 


32 


122.59694 


42.2739 


0.064 


-22.719 


2.40 


1.013 


3.7 


4.56 


5.02 


2.4 


2.9 


0.10 


-0.321 


-1.31 


1.298 


2.571 


0.9 


0.63 


0.21 


2 


33 


184.62019 


42.4610 


0.073 


-22.029 


1.93 


0.925 


-0.4 


9.19 


9.73 


5.1 


4.6 


-0.05 


-0.054 


-1.23 


-1.015 


2.284 


0.9 


0.15 


0.17 


2 


34 


254.08789 


39.2752 


0.062 


-22.595 


2.42 


0.976 


3.2 


294.65 


50.08 


10.7 


10.1 


-1.20 


-0.700 


-1.01 


-0.038 


2.422 


0.6 


0.48 


0.31 


2 


35 


129.46447 


26.5988 


0.088 


-22.705 


1.88 


0.903 


2.2 


6.32 


8.82 


6.2 


6.3 


-0.39 


-0.069 


-0.99 


-2.875 


2.413 


0.5 


0.92 


0.09 


2 


36 


134.49123 


30.3439 


0.086 


-22.817 


2.26 


0.949 


1.2 


15.11 


14.76 


7.4 


6.5 


-0.10 


-0.147 


-0.57 


-1.422 


2.439 


0.9 


0.11 


0.41 


2 


37 


133.51888 


29.0535 


0.084 


-22.938 


2.50 


1.002 


4.6 


8.06 


13.29 


2.9 


3.6 


1.72 


-0.099 


-0.49 


0.964 


2.463 


0.9 


0.35 


0.36 


2 


38 


237.97868 


27.8640 


0.082 


-21.860 


1.81 


0.979 


2.8 


8.80 


13.18 


2.9 


4.1 


-0.01 


-0.327 


0.29 


0.421 


1.996 


0.7 


0.56 


0.29 


2 


39 


136.68638 


3.6002 


0.072 


-22.123 


1.90 


0.893 


-1.5 


12.94 


21.63 


10.9 


10.2 


0.97 


0.146 


-0.51 


-1.928 


2.209 


0.5 


-0.07 


0.05 


2 


40 


144.40096 


7.9180 


0.093 


-21.778 


1.75 


0.925 


-1.3 


2.70 


4.25 


4.2 


4.4 


0.91 


-0.050 


-1.40 


-0.824 


2.270 


0.5 


0.20 


0.45 


2 


41 


141.65320 


3.4606 


0.088 


-22.235 


2.15 


0.954 


3.3 


4.45 


7.72 


4.8 


5.9 


-0.06 


-0.075 


-1.24 


-0.427 


2.422 


0.9 


0.21 


0.19 


2 


42 


178.15655 


3.4727 


0.081 


-22.771 


2.22 


0.964 


3.0 


5.71 


9.31 


1.6 


2.9 


0.41 


-0.625 


0.03 


-0.665 


2.417 


0.8 


0.54 


0.50 


2 


43 


7.36845 


-0.2126 


0.060 


-22.322 


1.85 


0.981 


5.3 


4.82 


4.83 


2.7 


2.5 


0.47 


0.019 


-1.35 


0.288 


2.390 


0.8 


0.41 


0.32 


2 


44 


216.19754 


2.6644 


0.053 


-22.400 


1.95 


0.961 


3.4 


7.04 


6.69 


4.5 


4.5 


1.04 


-0.009 


-0.90 


-0.344 


2.401 


0.8 


0.34 


0.25 


2 


45 


223.93717 


-0.3062 


0.083 


-22.295 


1.92 


0.905 


1.7 


1.09 


1.49 


2.0 


2.4 


0.91 


-0.369 


-0.81 


-1.818 


2.440 


0.8 


0.24 


0.20 


2 


46 


350.73485 


-10.0456 


0.084 


-22.368 


2.24 


0.989 


5.6 


4.94 


10.10 


5.6 


7.8 


1.13 


0.015 


-1.10 


0.477 


2.444 


0.9 


0.32 


0.36 


2 



